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FOREWORD 


The  principal  aim  of  Contract  F-19628-78-C-0032  is  to  elucidate 


the  influence  of  water  as  an  atmospheric  constituent  on  the  over-all 
behavior  of  the  atmosphere.  The  processes  involving  water  which  will 
particularly  concern  us  are  evaporation  from  the  earth's  surface  with  its 
accompanying  cooling  of  the  surface,  condensation  or  evaporation  within 
the  atmosphere  with  its  accompanying  warming  or  cooling,  precipitation, 
reflection  of  short-wave  radiation  by  clouds,  and  absorption  and  emission 
of  long-wave  radiation  by  clouds  and  water  vapor.  We  are  interested  in 
how  the  atmospheric  wind  and  temperature  fields  differ  from  what  they 
would  be  if  water  were  not  present,  and  in  how  the  various  processes 
which  we  have  mentioned  contribute  to  the  differences.  We  are  also 
interested  in  such  matters  as  the  contributions  of  these  processes  to 
the  changes  in  the  wind  and  temperature  fields  which  would  accompany 


a  change  in  the  intensity  of  the  solar  radiation. 

Our  specific  goal  has  been  the  construction  of  a  low-order  numerical 

model  of  the  general  atmospheric  circulation  which  incorporates  all  of 

the  processes  involving  water  which  we  have  mentioned.  By  a  low-order 

model  we  mean  one  defined  by  a  minimal  number  of  prognostic  ordinary 

differential  equations  which  can  yield  a  qualitatively  valid  description 

of  the  processes  involved.  We  do  not  expect  a  low-order  model  to  give  ~or  ~~ 

quantitatively  correct  results,  and  we  certainly  do  not  expect  it  to  ^ 

d  pi 

produce  good  day-to-day  weather  forecasts.  f J 

The  most  obvious  advantage  of  a  low-order  model  is  the  speed  with 

which  numerical  integrations  can  be  performed.  This  allows  us  to  examine  °13/ 
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many  more  cases  than  could  economically  be  treated  with  a  large  model. 
Perhaps  equally  important  is  the  fact  that  the  paucity  of  dependent 
variables  greatly  facilitates  the  analysis  of  the  results. 

The  principal  section  of  this  report,  entitled  "The  feasibility  of 
a  low-order  model  of  a  moist  general  circulation",  is  devoted  to  a 
description  of  the  low-order  moist  model  which  we  have  formulated.  We 
also  describe  some  preliminary  numerical  experiments  in  which  the  solar 
heating  is  taken  to  be  horizontally  uniform.  We  intend  to  submit  this 
section  to  a  journal  for  publication,  af‘  r  we  have  performed  and  analyzed 
a  few  numerical  experiments  with  horizontally  variable  solar  heating. 

Under  a  previous  contract  (F  19628-77-C-0026)  we  introduced  the 
concept  of  moist  available  energy  (MAE).  This  is  identical  to  available 
potential  energy  in  a  dry  atmosphere,  but  in  a  moist  atmosphere  it  is 
defined  in  such  a  manner  that  both  dry-adiabatic  and  moist-adiabatic 
processes  preserve  the  sum  of  MAE  and  kinetic  energy.  MAE  may  be  generated 
and  destroyed  by  diabatic  heating  and  cooling.  It  may  also  be  generated 
by  evaporation  from  the  earth's  surface  and  destroyed  by  precipitation, 
but  it  is  not  altered  by  condensation  and  evaporation  within  the  atmosphere. 
We  devised  a  graphical  procedure  for  evaluating  the  MAE  present  in  a  given 
atmospheric  state. 

One  of  the  first  tasks  under  the  present  contract  was  to  devise  a 
numerical  procedure  for  evaluating  MAE.  We  felt  at  that  time  that  this 
task  should  precede  construction  of  our  low-order  model,  and  that  the 
low-order  model,  when  formulated,  should  treat  MAE  in  a  qualitatively 
correct  manner.  As  our  low-order  model  stands  at  present,  we  have  not 


identified  any  quantity  which  possesses  the  properties  of  MAE. 
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Having  developed  a  rapid  numerical  procedure  for  evaluating  MAE, 
we  felt  that  it  would  be  desirable  to  apply  it  to  a  sequence  of  real 
weather  situations.  In  order  to  perform  the  evaluations,  we  had  to 
introduce  assumptions  regarding  the  amount  of  liquid  water  at  each  point 
in  the  atmosphere,  since  liquid  water  is  not  one  of  the  quantities  which 
is  observed  on  a  routine  basis.  We  soon  discovered  that  our  results  were 
highly  sensitive  to  the  manner  in  which  we  estimated  the  liquid  water 
content.  We  therefore  decided  to  abandon  our  attempts  to  study  the  day- 
to-day  variations  of  MAE,  at  least  until  such  time  as  suitable  methods 
of  estimating  liquid  water  content  became  available,  possibly  through  new 
satellite  observations. 

Another  of  our  aims  was  to  develop  a  low-order  moist  model  spanning 
many  scales  of  motion,  possibly  from  the  global  scale  to  the  cumulus 
scale,  in  contrast  to  the  model  described  in  this  report,  which  covers 
only  the  largest  scales.  It  became  evident  that  such  a  model  could  not 
easily  be  formulated  until  the  present  moist  model  was  in  working  order. 
As  a  result  we  have  not  progressed  with  this  task. 

As  described  in  the  main  section  of  this  report,  we  have  already 
uncovered  one  surprising  result,  namely  that  with  certain  intensities 
of  horizontally  uniform  solar  heating,  there  are  two  markedly  different 
equilibrium  distributions  of  temperature  and  moisture.  We  believe  that 
further  experiments  with  the  model  will  uncover  some  equally  interesting 


results. 
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THE  FEASIBILITY  OF  A  LOW-ORDER  MODEL 


OF  A  MOIST  GENERAL  CIRCULATION 


1.  Introduction 

The  past  two  decades  have  seen  the  ever-increasing  application  of 
hiqhly  simplified  systems  of  equations  to  various  problems  involving  the 
circulation  of  the  atmosphere  or  the  ocean  or  some  other  fluid  system. 

Such  systems  of  equations  have  come  to  be  known  as  "low-order  models." 

The  construction  of  a  low-order  model  of  the  atmosphere  usually  begins 
with  the  selection  of  a  set  of  n  prognostic  partial  differential  equations, 
and  some  auxiliary  diagnostic  equations,  which  approximate  the  equations 
governing  the  atmosphere;  the  n  dependent  variables  usually  include 
the  wind  components  and  the  temperature,  and  sometimes  the  water  content. 

The  vertically  continuous  atmosphere  is  first  approximated  by  an  array 
of  N  layers;  the  three-dimensional  field  of  each  dependent  variable 
is  replaced  by  fsl  two-dimensional  fields,  and  the  equations  are  modified 
accordingly.  Often  InJ  =  ;  sometimes  hi  =  \  .  Each  two-dimensional  field 

is  then  expanded  in  a  series  of  orthogonal  functions  -  generally  a  double 
Fourier  series  if  the  earth  is  approximated  by  an  infinite  plane  or 
channel,  or  a  series  of  spherical  harmonics  if  spherical  geometry  is 
retained  -  and  the  coefficients  of  the  orthogonal  functions  become  the 
dependent  variables  in  the  new  equations.  Finally  the  system  is  truncated 
by  omitting  reference  to  all  but  a  small  number  K  of  terms  in  each 
series.  The  resulting  xr'iN  prognostic  ordinary  differential  equations 
can  retain  some  of  the  essential  nonlinearity  of  the  advective  process 
even  when  K  -  3  • 


—  ■>  — 


A  similar  procedure  may  be  used  when  the  fluid  system  is  a  portion 
of  the  atmosphere,  all  or  a  portion  of  the  ocean,  or  a  laboratory  model 
of  the  atmosphere  or  the  ocean.  Variations  are  possible;  for  example, 
the  value  of  M  or  K  may  be  different  for  different  variables. 

Ordinarily,  the  advective  terms  in  a  low-order  model  will  preserve  total 
energy  if  those  in  the  original  equations  do  so  also,  provided  that  the 
total  energy  is  a  quadratic  function  of  the  dependent  variables. 

One  does  not  ordinarily  expect  quantitatively  correct  results  from 
a  low-order  model  when  the  truncation  of  the  series  is  severe.  It  is 
widely  assumed,  however,  perhaps  mainly  as  an  article  of  faith,  that  a 
carefully  constructed  low-order  mndei  can  give  results  which  are  qualitatively 
realistic,  at  least  over  certain  ranqes  of  the  adjustable  parameters. 

The  advantages  of  low-order  models  are  that  in  some  instances  analytic- 
solutions  may  be  found,  while,  when  numerical  methods  must  be  used,  the 
computations  needed  to  simulate  a  given  span  of  real  time  may  be  fewer 
by  several  orders  of  magnitude  than  those  required  when  a  large  global 
circulation  model  or  something  equally  detailed  is  used.  Also,  the 
paucity  of  variables  can  facilitate  the  subsequent  interpretation  of  the 
numerical  results.  Thus,  with  a  model  consisting  of  only  eight  prognostic- 
ordinary  differential  equations,  we  were  able,  using  analytic  procedures, 
to  reproduce  the  qualitative  arrangement  in  parameter  space  of  the  transi¬ 
tion  between  the  Hadley  and  Rossby  regimes,  as  observed  in  laboratory 
experiments  (Lorenz,  1962).  When  we  increased  the  number  of  equations 
to  fourteen,  we  were  able,  using  numerical  procedures,  to  reproduce  the 
transition  from  steady  Rossby  flow  to  vacillation,  with  relatively  little 
computational  effort  (Lorenz,  1963) . 
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Many  of  the  recent  large  global  circulation  models  treat  the  atmosphere 
as  a  mixture  of  dry  air,  water  vapor,  liquid  water,  and  perhaps  ice. 
Low-order  models  of  small-scale  systems  also  have  frequently  been  moist 
models.  Low-order  qeneral-eirculation  models,  however,  have  generally 
treated  the  atmosphere  as  an  ideal  gas.  Reasons  for  this  state  of  affairs 
are  not  hard  to  identify.  First,  although  simple  products,  such  as  those 
appearing  in  advective  terms,  assume  a  reasonably  simple  form  when 
transformed  into  Fourier  space,  some  of  the  nonlinear  functions  arising 
in  moist-atmosphere  thermodynamics,  such  as  those  relating  saturation 
mixing  ratio  to  temperature  and  pressure,  transform  into  very  complicated 
expressions . 

With  large  models  this  difficulty  may  be  circumvented  by  transforming 
the  variables  from  Fourier  space  to  grid-point  space,  performing  the  thermo¬ 
dynamic  computations  at  each  grid  point,  and  then  transforming  back  to 
Fourier  space,  at  every  time  step,  thanks  to  the  Fast  Fourier  Transform. 

With  low-order  models  this  procedure  introduces  another  difficulty.  The 
mixing  ratio,  which  is  commonly  used  as  a  measure  of  water-vapor  content, 
is  always  positive,  but  on  a  global  scale  it  may  vary  by  one  or  even  two 
orders  of  magnitude  in  the  horizontal  direction.  A  highly  truncated 
Fourier  series  which  best  approximates  the  mixing  ratio  over  the  bulk 
of  the  atmosphere  is  likely  to  indicate  negative  values  at  some  points 
and  value;-,  representing  suporsuturation  at  others. 

The  purpose  of  this  work  is  to  demonstrate  that  the  formulation  of 
a  low-order  model  of  a  moist  general  circulation  is  nevertheless  feasible. 

We  shal  l  construct  a  model  in  which  the  basic  variables  are  Fourier 
coefficients,  and  we  shall  transform  the  variables  to  and  from  grid- 
joint  space.  Wo  shall  include  the  thermodynamic  effects  of  phase  changes 
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of  water,  and  also  the  effects  of  water  vapor  and  liquid  water  on  short¬ 
wave  and  long-wave  radiation,  which  we  shall  formulate  in  a  simplified 
manner  compatible  with  the  low  vertical  resolution  of  the  model.  We  shall 
not  include  the  possible  presence  of  ice.  We  shall  allow  liquid  water 
to  remain  for  a  while  in  the  form  of  clouds  before  returning  to  the 
earth  as  rain. 

Since  the  atmosphere  receives  much  of  its  energy  in  the  form  of  latent 
heat  through  evaporation  from  the  ocean,  we  shall  find  it  convenient  to 
let  the  earth's  surface  consist  entirely  of  ocean.  The  simplest  arrange¬ 
ment,  in  the  sense  of  entailing  the  fewest  dependent  variables,  would 
be  to  let  the  ocean  have  zero  heat  capacity,  and  to  determine  its  temp¬ 
erature  diagnostically  by  requiring  that  it  neither  gain  nor  lose  energy. 
Such  an  ocean  might  more  properly  be  thought  of  as  wet  land.  It  is  more 
realistic,  and  just  as  simple  computationally,  to  give  the  ocean  a  finite 
heat  capacity,  in  which  case  the  ocean  temperature  becomes  another  depen¬ 
dent  variable.  We  shall  not  consider  an  ocean  of  infinite  heat  capacity, 
which  would  require  prespecifying  the  ocean  temperature. 

To  avoid  the  particular  difficulties  which  we  have  noted,  we  shall 
make  a  few  departures  from  the  most  commonly  used  procedures.  First,  we 
shall  use  a  quantity  representing  total  water  content  rather  than  water- 
vapor  content  as  a  basic  variable.  We  shall  use  an  auxiliary  diagnostic 
equation  to  specify  how  much  of  the  total  water  is  vapor  and  how  much  is 
liquid;  if  this  equation  is  judiciously  formulated  the  problem  of  super¬ 
saturation  will  never  arise.  In  formulating  the  equation  we  shall  recog¬ 
nize  that  each  grid  point  actually  represents  a  large  area,  and  that  a 
portion  of  the  area  may  be  subsaturated  while  another  portion  contains 


clouds . 
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Next,  our  variable  representing  total  water  content  will  be  the 
"total  dew  point"  rather  than  the  total-water  mixing  ratio.  By  the  total 
dew  point  we  mean  the  value  which  the  ordinary  dew  point  would  acquire 
if  all  the  liquid  water  were  converted  to  vapor,  at  constant  pressure. 

Like  the  temperature,  the  total  dew  point  should  vary  by  a  factor  of 
less  than  two  in  the  horizontal  direction,  and  even  a  highly  truncated 
Fourier  series  representing  it  should  not  produce  negative  values. 

In  the  following  sections  we  shall  formulate  first  the  equations  of 
a  continuous  atmosphere,  then  the  equations  of  a  layered  atmosphere,  and 
finally  the  equations  of  a  low-order  model.  Since  we  shall  be  concerned 

mainly  with  feasibility,  we  shall  at  times  suggest  alternative  formulations, 
and  out  numerical  values,  when  not  dictated  by  the  physics,  will  simply 

be  suggested  values.  Wo  shall  present  the  results  of  a  few  numerical 

integrations,  which  will  constitute  an  initial  test  of  the  model. 


2.  The  continuous  equations 

Our  independent  variables  will  be  time  't  ,  pressure  p  ,  and 
horizontal  coordinates  X  and  .  Our  basic  dependent  variables 

will  be  heiqht  7  ,  horizontal  velocity  derivable  from  a  stream  function 

4'  and  a  velocity  potential  X  ,  individual  pressure  change  U!  , 
total  dew  point  W  ,  air  temperature  "T  ,  and,  in  our  principal 
version,  sea-surface  temperature  S>  -  From  these  we  shall  derive  the 
auxiliary  variables  total-water  mixing  ratio  mt  ,  water-vapor  mixing 
ratio  a r  ,  saturation  mixing  ratio  M  corresponding  to  air  temperature 

and  pressure,  and  saturation  mixing  ratio  S  corresponding  to  sea- 
surface  temperature  and  sea-level  pressure.  We  shall  let  the  atmosphere 
occupy  a  channel  of  infinite  west-east  extent,  bounded  laterally  by  fric- 

,  and  below  by  the  surface  *  p 


tionless  walls  at  -  0 


and 


whose  height  will  be  allowed  to  vary.  Suggested  values  for  D  and 
Pg  are  2000  km,  if  the  channel  is  supposed  to  simulate  the  extra- 
tropical  latitudes  of  one  hemisphere,  and  1000  mb. 

Our  basic  diagnostic  equations  will  be  the  hydrostatic,  qcostrophic, 
and  continuity  equations 


dz/sp  =  -8T/(}q  , 


(1) 


f.t  =  3- 2  ■ 


(2) 


-V  =  o 


(3) 


where  R>.  is  the  gas  constant  for  air,  cp  is  the  acceleration  of  gravity 
and  f o  is  the  constant  average  value  of  the  variable  Coriolis  parameter 
•p  .  We  shall  let  CO  vanish  at  the  top  and  bottom  of  our  model 
atmosphere . 

Our  prognostic  equations  will  be  the  vorticity  equation 


(4) 


the  equation  for  continuity  of  water  content 


M  *r  j  J  t  -  G-  J 


(5) 


and  the  first  law  of  thermodynamics  applied  separately  to  the  atmosphere 


and  the  ocean, 
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-  RTw/p  +  H  ,  (6) 

J  (cS)  /  J-t  =.  E  (7) 

where  L_  is  the  latent  heat  of  condensation  of  water  vapor,  assumed 

constant,  Cp  is  the  specific  heat  of  air  at  constant  pressure,  C  is 

the  specific  heat  of  liquid  water,  and  X  denotes  a  Jacobian  with  respect 

x  _ 

to  and  .  Here  V  F  J  G-,  H  ,  and  t  are  source  and  sink 

terms  for  vorticity,  atmospheric  water  content,  and  atmospheric  and  oceanic 
enthalpy;  jp"  includes  friction,  G"  includes  evaporation  and  precipitation, 
and  H  and  E  include  diabatic  heating  and  the  effects  of  evaporation 
from  the  ocean  surface  (but  not  the  effects  of  condensation  within  the 
atmosphere) .  By  writing  (6)  as  a  prognostic  equation  for  specific  enthalpy 
crT  +  L  ^  ,  we  automatically  include  the  thermodynamic  effects 

of  water.  The  omission  from  (4)  of  the  nonlinear  terms  containing  'X- 
or  CO  is  consistent  with  the  geostrophic  approximation.  Suitable 

numerical  values  for  the  constants  arc  "f0  =  10  ^ 9.8  ms~X 

Cp=  1000  ms  K  ,  C  =  4185  m  s  K  ,  P.  =  (‘i./'7)C  »  and  L.  =  2.5  x  10  m  s 


Wo  shall  relate  the  mixing  ratios  to  the  temperatures  by  the  approxi¬ 
mate  formulas 


A*'"  -  C*  W7  P  i 


(8) 


(9) 


(10) 
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I 

where  C.  and  ^  are  constants  for  which  appropriate  values  will 
presently  be  introduced.  To  close  the  system  we  need  an  equation  relating 
fj-  to  the  other  variables.  The  standard  assumption  is  that  there  is 
no  liquid  water  if  the  total  water  is  insufficient  for  saturation,  i.e., 
nt-  !+r  if  ,  and  no  supersaturation  if  the  total  water  is  suffi¬ 

cient,  i.e.,  N-  =  m  is  jur  >  w  .  Sometimes  condensation  is  assumed 
to  occur  when  the  relative  humidity  A  =.  at/ aa.  reaches  some  value  below 
saturation,  perhaps  80  per  cent.  As  we  have  noted,  we  prefer  a  formulation 
where  a  portion  of  the  region  represented  by  a  single  grid  point  may  be 
subsaturated,  while  another  portion  may  contain  clouds.  We  shall  choose 
a  formula  which  makes  the  liquid  water  content  /ur  ~  at-  small  when  the 
degree  of  subsaturation  M  -  N~  is  large,  and  vice  versa.  A  simple 
formula  of  this  sort,  which  makes  /wj w  -»  |  as  /*r  -»  q  ,  and  *r / m  -*  \ 
as  /*r  -)  oc  »  is 


(/^r  -  tr)  [  *  ~  AT  ^  -  yV" 


(11) 


where  Y  is  a  small  constant.  Eqs.  (1) -  (11)  form  our  closed  system  in 
the  11  basic  and  auxiliary  dependent  variables. 

Choosing  V  a  '/h  makes  A  =  0.^7when  X  ’  A  =  O  .  9  O  when 

fijj-  _  ^  ,  and  r\_  s  0  *1  5  when  M  •  We  also  note  that  choosing 

Y  as  0  reduces  (11)  to  the  standard  assumption. 

To  arrive  at  formulas  (8)- (10)  we  write  the  Clausius-Clapeyron  equation 
in  the  approximate  form 


t)/JT  =  t  L  a.*  1*0/  (  R  T*t  ) 


(12) 


A 
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where  J2S  M  is  the  saturation  vapor  pressure  of  water  at  temperature  T  , 
t  =0.6^^  is  the  ratio  of  the  molecular  weights  of  water  and  air,  and  T 
is  a  constant  temperature  typical  of  the  atmosphere,  say  273  K.  This 

_*  -r 

modification  from  the  more  exact  form,  in  which  I  instead  of  I  I 
would  appear  in  the  denominator,  and  L  would  vary  with  T  ,  makes 
proportional  to  T"^*  »  where  /a  =  a/(dT")  is  a  constant 

whose  appropriate  value  is  about  20.0.  Eq.  (8)  then  follows  from  the 
approximate  relation  AA  r  &s ( ”0  /  P  ;  the  constant  c*  is  to  be 
chosen  so  that  M  =  0.0038  when  f  =  273  K  and  p  =  1000  mb.  Eqs. 

(9)  and  (10)  are  arrived  at  similarly. 

From  (1)  and  (2)  we  obtain  the  thermal  wind  equation 


3^/  i?  =  -  RT/(f.p) 


(13) 


after  which  further  reference  to  Z  is  superfluous.  Our  closed  system 
is  now  (3)- (11)  and  (13). 

We  would  prefer  a  prognostic  equation  for  w  to  one  for  mJ~ 

From  (5)  and  (8)  we  find  that 


/(W^t  -  VV/toy'p  +VWG-//W-  t 


(14) 


where  )/  »  '/r  Likewise,  we  would  prefer  a  prognostic  equation  for  T 
to  one  for  C-^T  +  L  W"  .  We  find  from  (11)  that 


1 


-13- 


where 


NT„  =  Ar  (M  '  ^/[a*t(m-  ^  +  M  (w  -  a r)~)  # 
Ar^  -  /v-  (/*r-Ar^y/  £>*«-(  A4-  )  -+  M  [  W  -  at)  ^ 


With  the  aid  of  (9)  and  (15),  (6)  becomes 


1  Cr+c*VT/n  -  (xcrtvct)rw/  p  -b(H-L^-G-) 

where  R/c,  and  C.  q_  z  /*  L  M  /  T  . 

Finally,  we  need  to  rewrite  (17)  and  (18)  with  local  time 
We  obtain 

2w/n  --  -tM  -  VX'VW  -  w 

+  v  +  v  w  (r  [  a*7" 

JT/It  =  -J^T)  -  Vl'VT 

+  [(K^VCi)/h«t)]Tu/p  +(H  -L  <>)/( 


(15) 


(16) 


(17) 


(18) 


derivatives. 


(19) 


Cr+Cf).  (20) 
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To  interpret  the  coefficient  of  Toj/|p>  in  (20)  we  may  picture  an 
adiabatic  chart,  with  pressure  and  temperature  as  coordinates,  drawn  so 
that  the  slope  of  the  dry  adiabats  is  X  .  The  slope  of  the  lines  of 
constant  saturation  mixing  ratio  is  then  V  ,  while  the  ratio 

,  which  lies  between  X.  and  \)  ,  is  the 

effective  slope  of  the  moist  adiabats.  It  would  equal  the  standard  slope 
of  the  moist  adiabats  if  *  l  ,  and  that  of  the  dry  adiabats  if 

-o  -  The  factor  a J"M  in  c  9-  is  defined  by  (17),  and  actually 

;  |  if  A.  =■  |  and  =  0  if  A  -  O  ,  and  0  <  AT^  |  otherwise. 

Eq.  (20)  thus  recognizes  that  under  adiabatic  lifting  f\,  increases, 
so  that  there  is  somewhat  less  condensation,  and  the  rate  of  cooling  is 
somewhat  closer  to  the  dry-adiabatic,  then  if  T L  remained  fixed  at 
100  per  cent. 

3.  The  layered  model 

For  our  layered  model  we  shall  choose  a  simple  form  of  the  two-layer 
model,  in  which  the  basic  dependent  variables  include  ip  at  each  of  two 

levels,  but  T  at  only  one  level.  Such  models  have  been  called 

"2-1/2-dimensional",  since  the  wind  field  is  3-dimensional  while  the  temp¬ 
erature  field  is  effectively  2-dimensional.  It  is  consistent  with  this 
formulation  to  lot  the  basic  variables  include  w  at  only  one  level. 

We  shall  let  the  atmosphere  consist  of  an  active  troposphere  extending 
upward  from  1000  to  200  mb,  and  an  isothermal  dry  stratosphere  extending 
upward  from  200  mb.  The  stratosphere  will  play  no  role  except  to  exert 
the  200-mb  pressure  at  its  base.  We  shall  let  the  ocean  consist  of  an 
isothermal  mixed  layer  extending  downward  to  8000  mb  (about  70  m) ,  and  a 
deeper  ocean  which  will  play  no  role.  We  shall  use  subscripts  0,1, 2, 3  and 
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4  to  denote  the  value  of  an  atmospheric  quantity  at  200,  400,  600,  800, 
and  1000  mb  respectively,  and  subscripts  5  and  6  to  denote  an  oceanic 
quantity  at  4500  and  8000  mb.  We  shall  sometimes  choose  other  numerical 
values  for  and  |°fc 

Our  basic  dependent  variables  will  then  be  'fj  ,  , 

,  W  ,  Ta  ,  and  S  s  .  We  must,  however,  recognize  the  existence 

of  the  variables  at  other  levels,  and,  in  formulating  expressions  for 

radiative  heating,  we  must  let  Vv'j.  and  determine  W  and  T  at 

other  levels  in  a  reasonable  manner.  We  shall  first  require  that  in  each 

vertical  column,  T  be  proportional  to  ,  where  X  is  a  prespecified 

constant.  Noting  that  X  a  O  or  X  ~  ^  -  ^~/l  would  imply  respectively 

an  isothermal  or  dry-adiabatic  lapse  rate,  we  find  that  a  reasonable  value 

of  \  would  be  about  0.175.  We  next  specify  that  the  relative  humidity 

/I  be  constant  within  a  column.  It  then  follows  that  hr j  m.  ,  and 

hence  W/  T  is  constant,  so  that  w  is  also  proportional  to  Ir*  , 

(Ay m-0 

while  ,  at  ,  and  M  are  proportional  to  p  / 

From  (3)  and  (13)  we  find,  since  U)0  =  U)^  ■=.  O  ,  that 

V1*  ,  (21 ) 


vV,  ■  ,  ,22) 

=  Vi  t  ,  _  1231 

where  -i °.V  *  and  Vi  s  RA  p/tt.pO  .  In  view  of  (21)  and  (22), 

the  vorticity  equation  (4)  yields 
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+  i-7VRv  7 

-'  -*  0 -^^/ap  +  j  . 

Letting  %*['*<+%)/?■  ,  we  find  from  (23)-  (25)  that 


--  -TK/rn)  -  Tl^vV^f)  +  }R«>s/(hV^ 
+  vx(f,-f5)/k  , 

while 


-^jKyxo/M  +v(f,+f1)/i  (27) 


Because  w  and  "T  are  proportional  to  p  , 
and  Cv>  in  (19)  and  (20)  may  be  replaced  by 

y  T  ^  l  f*  ,  and,  since  we  may  assume  that  Xa 

(19)  and  (20)  yield 


the  terms 

W  »>  j  p 


to  aw/ap 

and 


z  =  Oj 


+  (V->)  +  v  ^  ^ 

2  1>./zt  -*■  [(xc^yc^O/U^c^^T^/p^ 

+  [h%  -  l  c**')/ ■*■ 

*  (29) 
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while  (7)  becomes 

5Ss/>t=  Es/c 


The  diagnostic  equations  (8)- (12),  with  suitable  subscripts  added,  complete 
the  system. 

We  may  now  eliminate  the  time  derivatives  from  (26)  and  (29) ,  and 
obtain  the  to  -equation 

-  -  V>,  +  f) 

(31) 

a 

where  V  is  the  inverse  of  V  .  With  Ud  given  by  the  diagnostic 

equation  (31)  ,  we  end  up  with  a  system  of  four  prognostic  equations 
(27)- (30)  governing  the  dependent  variables  v|/.5  T -j.  ^  and 


4.  Sources  and  sinks 

Although  the  new  dependent  variables  are  ostensibly  the  values  of 
quantities  at  specific  levels,  they  are  supposed  to  represent  entire 
layers,  and  their  changes  due  to  external  influences  should  represent  the 
effect  of  the  external  influences  upon  entire  layers.  For  example,  we 
wish  T*  to  be  affected  by  the  total  absorption  and  emission  of  radia¬ 
tion  by  the  atmosphere,  rather  than  just  the  convergence  of  radiative  ■ 

j 

flux  at  level  2,  which  would  be  difficult  to  estimate.  i 


£ 
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J-  _  i.  _ 

The  quantities  ^  Ft  ,  V  F$  ,  tr  j  ,  H  ^  ,  and  t  ^ 


represent 


sources  of  vorticity  of  momentum,  water,  and  enthalpy  per  unit  mass. 

On  the  other  hand,  such  quantities  as  radiative  flux  and  exchanges  across 
the  ocean-atmosphere  interface  are  generally  expressed  in  amounts  per 

unit  horizontal  area.  Accordingly,  we  shall  denote  the  respective  sources 

a  /  *  »  ,  1  c  1 

per  unit  area  by  V  r,  ,  V  ^  ,  “j  ,  and  C5  .  We  find 


t _ .  —  l  p  i 

F  ,=  fir  %  F, 


Fa  -  A  r'1  %  Fj  , 


and,  in  view  of  the  assumed  vertical  variations  of  and  T  , 


>/*-'/  >r\"‘  r  ' 

g-  j  -  >/<  fy  In  -  r*  I  , 

H  J  (pm  -  P*  3  9-  H  x  # 


while 


I  i 


Es  -  (Pi-N)  j  Es 


For  such  a  crude  model  a  nonlinear  formulation  of  friction  seems 
unwarranted,  and  we  shall  let 
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F,'  =  -  h'  (+.- %)  , 

F,‘  *  -  k<f\,  +  *>’(  't-  , 


where  (3  ^  ^  .  An  appropriate  value  for  ta.  would 

be  Cx,  />4  U„  ,  where  Cj>  is  a  surface  drag  coefficient,  say 

2.0  x  10  ,  and  and  U^.  are  a  typical  density  and  wind  speed 

at  level  4,  say  1.2  kg  m  3  and  10  m  s  Probably  all  that  can  be  said 

about  K  is  that  is  should  be  smaller  than  k  ,  and  we  shall  let 

k'  =  k/ <4 

The  processes  entering  Gr  x  are  evaporation  and  precipitation.  We 
shall  let 


Crx  *  W  (  sr  *  “  k  *  A/0  , 


where  =.  ,  and  k  is  the  same  as  in  (38)  .  In 

parameterizing  precipitation  we  have  assumed  a  "half  life"  for  clouds. 

We  observe  that  when  (39)  is  substituted  into  (34),  the  coefficient  of 

2^.-  ,  with  the  suggested  numerical  values,  becomes  about  (5  days)  3. 

If  the  half  life  for  clouds  is  about  16  hours,  the  e-folding  time  is  about 
.  n 

1  day,  and  ft  5  5  R 

Exchanges  of  sensible  and  latent  heat  across  the  atmosphere -ocean 
interface  may  be  handled  similarly.  We  find  that 


-*■'  mm  rm upjwswwip 


Hjl  -  ^  (S,  "  +•  k  L  (  S5  “  ■+  ^  3k  J  *4C 

^  3-KC,(SrTH)  ^L(£f-^0  ,  (41 

>  »  I 

where  =  *^a  (|“h  /  Pa)  ,  and  Kj_  and  K5  denote  the  effects 

of  radiation. 

We  shall  assume  that  all  short-wave  radiation  which  is  not  reflected 
by  clouds  is  absorbed  by  the  ocean.  Short-wave  radiation  then  does  not 
enter  R^  ,  and  its  contribution  to  R  5  is  G?e(l  -  O')  ,  where  Q  0 

is  the  intensity  of  the  solar  beam  and  is  the  albedo.  We  shall  let 

clouds  be  perfect  reflectors  of  short-wave  radiation,  whence  the  albedo 

equals  the  fraction  of  the  sky  covered  by  clouds.  This  fraction  may  be 

H 

parameterized  in  a  number  of  ways;  the  formula  <X  -  JX-  seems  to 
yield  reasonable  values. 

We  shall  let  the  ocean  surface  be  a  black-body  absorber  and  emitter 
of  lonq-wave  radiation.  Absorption  and  emission  by  the  atmosphere  is  more 
complicated.  We  need  to  know  the  atmosphere's  fractional  absorptivity 
and  emissivity,  and  the  effective  temperatures  at  which  it  radiates  upward 
and  downward.  It  would  be  pointless  to  introduce  all  of  the  intracacies 
of  the  .appropriate  radiation  formulas  in  a  model  which  is  so  crude  in 
other  respects.  Instead  we  shall  use  a  modification  of  Simpson's  method. 

In  his  computations  of  the  atmosphere's  radiation  balance,  Simpson 
(1928)  treated  the  atmosphere  as  beinq  completely  transparent  in  the 


8.9-11  micron  band,  so  that  the  earth's  surface  rather  than  the  atmosphere 
would  radiate  to  space  in  these  wave  lenqths.  He  assumed  that  a  layer  of 


air  containing  0.3  nun  of  prccipi table  watei ,  i.a.,  0.3  kg  of  water  vapor 

2 

per  m  of  cross  section,  would  be  completely  opaque  below  7  and  above 
14  microns,  so  that  the  uppermost  such  layer,  generally  in  the  upper 
troposphere,  would  radiate  to  space  in  these  wave  lengths.  In  the  inter¬ 
mediate  bands  the  atmosphere  would  radiate  to  spaoqaL  an  intermediate 
temperature . 

We  shall  simplify  Simpson's  treatment  by  assuming  as  a  first 
approximation  that  the  cloud-free  fraction  of  the  atmosphere  radiates 

t 

upward  and  downward  with  a  certain  fraction  (X.  ,  say  0.7,  of  the  intensity 

of  black-body  radiation,  at  temperatures  and  respectively. 

/ 

We  shall  let  the  cloud-free  atmosphere  absorb  the  fraction  (X  of  the 

i 

radiation  from  the  earth's  surface,  while  the  fraction  I  -  passes 

I 

through.  We  shall  neglect  any  variations  of  <X  which  ought  to  occur 
because  the  wave  length  of  maximum  absorption  and  emission  shifts  with 

temperature.  We  shall  let  the  clouds  radiate  upward  and  downward  as 

i  ^  l 

black  bodies,  also  at  temperatures  T0  and  1 g  .  These  temperatures 

1  i 

will  occur  at  pressures  po  and  p^  ,  which  will  be  the  levels 
above  which  and  below  which,  respectively,  the  amount  of  water  vapor 


V  j  3-  ,  where  V  =  0.3  kg 


Letting  rr  be  the  value  which  At" ^  would  possess  if  the  total 


water  vapor  in  a  column  were  V 
k*  -  I 

to  '  ,  that 


we  find,  since  at  is  proportioned 


(P.'f  -  r-*  r.^\  , 

•  p/'"  -  i  -  p.'r)  . 


1  mipyj* 't'wrww*-* 
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With  the  numerical  values  we  have  been  using,  at  is  about  3  x  10  kg  m 
Eqs .  (42)  and  (43)  appear  satisfactory  when  *-T is  large,  but  they 

/  i 

become  unreasonable  when  bT ^  is  very  small,  since  K  ?  P*  when 

i  i 

<  W  ,  and  and  move  right  out  of  the  troposphere  when 

*1*1  Ar  /2.  .  it  would  be  more  reasonable  to  have  and  ps 

i 

approach  an  intermediate  value  Pi  as  -i  O  -  Formulas  which 

fulfill  this  condition,  but  differ  only  slightly  from  (42)  and  (43) 
when  AT^  is  large,  are 


1 


') 


J 


(44) 


V, 


+  fjr 


(45) 


With  our  chosen  numerical  values  the  exponent  is  3.5  . 

In  addition,  in  the  limit  as  aT^  o  there  should  be  no  absorption 
or  emission  by  water  vapor.  Wc  shall  fulfill  this  condition  by  modifying 
the  factor  O.  by  the  additional  factor  A Tj.  /  +  At  ")  .  We  find, 

then,  that 


R3  --  +  ('-*)  «■'  +  <r  (Sj-  Tv'4  -T.'^ 


(46) 


R  -  -  crS5H  +  +  (•  '*-)  +  ^  *  Qo  ( 1 "  ^  # 


(47) 
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l  .  t  h  1  ,  > 

where  T.  -  (  y»0/pj}  T#  =  T3l(  Py/|»a)  ^and  o~  =  5.67  x 

-8  -3-4 

10  kg  s  K  is  the  Stefan- Bo ltzman  constant. 


5.  The  truncated  equations 

To  convert  our  system  into  a  low-order  model,  we  shall  approximate 
each  dependent  variable  by  a  truncated  double  Fourier  series  in  JC.  and 
We  shall  use  the  same  truncation  which  wo  used  previously  in  a  study  of 
vacillation  (Lorenz,  1963) ,  to  which  the  reader  is  referred  for  details 
not  included  here. 

Recalling  that  the  domain  of  the  variables  is  an  infinite  strip 
with  boundaries  at  —  O  and  ^  ;  V  D  ,  we  let  y.e  -  */T>  and  =  */l) 
If  X  then  represents  any  one  of  the  variables  bO  ,  'f*  ,  W  ,  T  , 
or  S  ,  we  let 

6 


X  -  2.  x.<h 

i  »  o 


(48) 


where 


$o=  * 

§\  s  *  s,v> 
f  a  1  s'* 

$  4  =  COS  £ 

$  H  SIK 

=  3. 

CoS  3- 


(49) 


for  some  positive  constant  V\  .  In  our  numerical  work  we  shall  let  V\  -  . 


L 


We  observe  that  a  j 


satisfy  the  relations 


, 

r  -  <X;  J 


(50) 


(51) 


where  the  bar  denotes  an  average  over  the  whole  domain,  d  i.j,  is  the 
Kronecker  delta ,  and  <X0  -  0;  *■  ,  =  a  -  I  +  yi  ;  a  r  ‘■f  +  ^  j 

and  <X  t  ~  ^  .  Each  dependent-variable  field  thus  consists  of  an  overall 

average,  a  superposed  zonally  symmetric  portion  with  two  north-south 
modes,  and  a  single  superposed  wave  with  two  north-south  modes  and  an 
arbitrary  longitudinal  phase.  Wo  have  ordered  the  variables  so  that 
subscripts  divisible  by  three  refer  to  the  zonally  symmetric  field.  We 
include  the  effect  of  a  variable  Coriolis  parameter  F  by  letting 

hz  -  *•/*■*■ 

In  general  an  arbitrary  nonlinear  function  Z  of 
and  hence  such  a  function  of  the  dependent  variables,  cannot  be  expressed 
exactly  as  a  linear  combination  of  ,  and  further  approxima¬ 

tion  is  needed.  In  view  of  the  orthogonality  relation  (50),  the 
appropriate  value  of  is  2  •  If  2  is  a  simple  product, 

say  Z  -  XY  and  y  and  Y  have  the  form  of  (48),  we  readily  find  that 


tV 

f' 


2  =  X  X;Y. 

o  r-r 
>  *  0 

Z,  »  (XY^,  -  [(XY),v  -  (  Y  Y\.  ]  / 'r>- 

2 1  =  (yy)#1  ■  ~ 

Z  J  =  O^).*  -t[UY)„  t  (*Y)„  +(YY^t]/jl 

f,  .  (xr)„  +  (<Y \Jfx 

2s  +(<A,/^ 

Zt  =  ( X  Y)#t  -  (X.Y,  -t 

where  (*YV  stands  for  4  X^  Y  *  *  likewise,  if 


we  find  that 


and 


where  L<A  stands  for  and  t  =  *v,ra/(,5*5>). 

When  ?  represents  a  more  complicated  function  of  the  dependent 
variables,  such  as  one  of  the  high  powers  or  fractional  powers  appearing 
in  the  source-sink  terms  in  (27) -(31),  an  analytic  expression  for  ~2-  ^ 
would  be  prohibitively  complicated.  To  approximate  it  would  be 

preferable  to  evaluate  2.  and  ^  ^  at  each  of  a  well-chosen  set 
of  points  spanning  the  domain,  and  then  to  calculate  the  average  product 
of  ^  and  ^  K 

Accordingly,  wo  shall  introduce  a  grid  of  points  into  the  domain.  At 
each  time  step  we  shall  evaluate  each  dependent  variable  (except  ''Pa  ) 
at  each  point,  and  then  compute  the  source-sink  terms  in  (27)- (31)  separately 
at  each  point,  finally  multiplying  each  value  so  obtained  by  the  value  of 
(£>  .  at  the  same  point,  and  averaging.  Since  the  functions  are 

trigonometric,  this  procedure  is  nothing  more  than  transforming  from 
Fourier  space  to  grid-point  space  and  back  again. 

The  estimates  of  ~Z  $  ^  ,  where  2  is  a  source-sink  term,  will  be 

more  accurate  the  larger  the  number  of  points  chosen,  but  since  the  model 
as  a  whole  is  rather  crude  regardless  of  how  many  points  are  chosen,  there 


is  little  to  gain  by  choosinq  too  many.  We  shall  settle  for  a  grid  of 
16  points,  which  will  bo  located  at  the  intersections  of  the  lines 
#-o-  lf/*  J  3  ,  5>r/2j  1  ir/s  with  the  lines  v\  %0  =  ^  I  ^  j 

,  5  tr /^,  7^/4- 

To  solve  the  to  -equation  (31)  we  must  either  evaluate  a  product 

rt  3L 

implicitly  in  Fourier  space  or  invert  y  implicitly  in  grid-point 
space.  We  shall  choose  the  former  procedure.  We  first  evaluate  the 
coefficient  of  in  the  first  term  of  (31)  (call  it  0)  in  grid- 

f>oint  space,  and  then  transform  S  into  Fourier  space.  Since  HA 
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the  coefficients  of  ~  £  in  the  Fourier  transform  of  (31), 

with  (  ^  as  in  (52)  ,  form  a  set  of  six  algebraic  equations  to 

be  solved  for  - ,(0^)^  .  The  terms  (  in  the  Fourier 

transforms  of  (28)  may  then  be  evaluated,  using  (52). 

6.  Testing  the  model 

With  a  typical  low-order  dry  model  of  the  general  circulation  the 
determination  of  an  equilibrium  solution,  when  the  external  heating  is 
horizontally  uniform,  is  a  trivial  matter,  but,  with  our  moist  model  it 
is  not  at  all  obvious  what  such  a  solution  should  look  like.  The  motion 
should  vanish,  but  the  equilibrium  values  of  W  ,  T-,  ,  and  S  c 

are  solutions  of  rather  complicated  nonlinear  equations,  obtained  by 
equating  ,  Hj  ,  and  E  5  to  zero. 

As  an  initial  test  of  the  model  we  shall  determinethese  equilibrium 
values  for  various  values  of  Q  .  More  definitive  tests  with  horizontally 
varying  Q  have  yet  to  be  performed.  Assuming  that  a  computer  program 
to  integrate  the  equations  of  the  model  has  been  written  in  any  case,  the 
way  to  determine  the  equilibrium  solution  which  requires  the  least  addi¬ 
tional  programming,  although  not  the  least  computing,  is  to  run  the  model 
from  arbitrary  initial  conditions  until  a  steady  state  is  approached.  Since 
there  are  no  horizontal  variations,  transformations  from  Fourier  space 
to  grid-point  space  and  back  again  are  unnecessary,  and  the  computation 
may  be  performed  at  a  single  grid  point.  Because  of  the  large  heat 
capacity  of  the  ocean,  the  approach  to  equilibrium  may  require  many  years, 
ard  we  should  note  that  the  final  state,  as  opposed  to  the  process  of 
approaching  it,  is  independent  of  the  ocean’s  heat  capacity.  We  shall 


therefore  use  an  ocean  with  a  2-m  instead  of  a  70-m  mixed  layer,  whose 


heat  capacity  is  comparable  to  that  of  the  atmosphere. 

We  shall  express  the  solar  heating  in  terms  of  a  planetary 


temperature  ,  where  T 


We  do  not  include  the  albedo 


Ou  in  the  definition  of  T <q  <  since  cx,  is  not  prespecified. 

We  have  obtained  several  numerical  solutions  using  the  suggested 
numerical  values  of  the  constants,  with  values  of  "Tq  varying  from 
258  K  to  288  K.  We  note  that  Qo  -  340  Wm-  ^ ,  the  approximate  average 
value  for  the  earth,  when  T*  --  Vi*  K  .  We  have  found  it  convenient 
to  use  the  surface  values  and  instead  of  w  ^  and  Ta 

as  output,  although  Vj,  and  remain  the  working  variables.  In 

each  run  VA,  ,  ~T  4  ,  and  S5  all  equal  Tq  initially; 

would  b<'  the  equilibrium  value  for  £->5  if  the  greenhouse  effect  exactly 
canceled  effect  of  the  albedo.  We  have  used  1.5-hour  time  steps.  After 
two  years,  in  each  run,  except  one  where  Tq -  413  K  ,  when  five  years 
are  needed,  the  approach  to  equilibrium  becomes  so  nearly  exponential 
that  the  final  state  can  be  easily  extrapolated  by  hand  computation. 

Table  1  shows  the  results.  The  most  outstanding  feature  is  the 
sensitivity  of  the  equilibrium  solution  to  Tq  ;  when  T3  is  low, 

,  and  S>  5-  are  very  low,  and  when  Tq  is  high,  Wuj  , 
Ta  and  ^5  arc  vory  high.  The  change  as  T"q  changes  from 
268  K  to  273  K  is  especially  abrupt.  It  is  also  noteworthy  that 
is  below  V/i)  and  $>5  when  is  low,  and  above  W4  and  Sf 

Yr  is  high,  while  Na/^  is  always  slightly  below 


when 
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Table  1.  Equilibrium  values  of  surface  total  dew  point  ,  surface 
temperature  ,  sea-surface  temperature  S,,.  and  relative  humidity  r, 
for  various  values  of  planetary  temperature  T{^,  as  determined  in  first 
set  of  numerical  runs. 


tq 

258 

263 

268 

273 

278 

283 

288 

W4 

224.1 

229.0 

234.5 

306.4 

315.9 

322.9 

329.1 

T4 

223.1 

227.6 

232.8 

312.4 

324.7 

334.2 

342.9 

S5 

225.4 

230.1 

235.7 

306.9 

316.3 

323.3 

329.4 

r 

0.84 

0.85 

0.85 

0.62 

0.53 

0.47 

0.42 

These  features,  once  observed,  are  readily  accounted  for.  Wo  note 
first  that  with  our  chosen  values  for  the  constants,  and  aT^ 

are  about  3.6  times  and  aT,  .We  thus  find  from  (39)  that  when 


& 


l 

3l 


vanishes. 


/WM  -  (  1  -  3-fe  k/k”)  ath  +-  (  3. 4.  fe  /k  )  s 

II 

It  follows  that  i f  we  had  chosen  k  =  2>,(>  k  ,  which  seems  just  as 
> 

reasonable  as  our  choice  of  5  K  *  and  S s  would  be  equal. 

With  fe.  =  5  k  ,  becomes  a  weiqhted  averaqe  of  AH#  and  S^-  , 

and,  since  <  /uTw  ^  S 5  >  ,  whence  S 5  r'  . 

Suppose  next  that  we  have  located  a  value  of  T for  which  -S5 
and  that  w e  then  increase  if  ,  Ty  ,  and 


w  <*ro 


to  increase  in  pro{X5rtion,  the  relative  humidity,  and  hence  the  cloud 


albedo,  would  remain  fixed.  The  latent  heat  transferred  from  the  ocean 
to  the  atmosphere  would  then  increase  as  ,  while  the  heat 

returned  to  the  ocean  through  radiation  or  sensible-heat  exchange  would 

4  -r- 

increase  only  as  I  ^  or  I  ^  Equilibrium  would  then  no  longer 

exist,  and  the  atmosphere  would  warm  up  while  the  ocean  cooled.  With 
equilibrium  reestablished,  ~T^  ^  5  g  •  Likewise,  if  Td  should 

decrease,  T 4  <•  Sy  . 


It  follows,  since  and  3^  vary  similarly,  that  the  relative 

humidity,  and  hence  the  cloud  cover  and  the  cooling  effect  of  its  albedo, 
must  decrease  as  %  increases.  On  the  other  hand,  the  greenhouse 
effect,  which  depends  ujK>n  the  total  water-vapor  content,  increases 
when  increases.  The  equilibrium  values  of  ,  "T^  ,  and  S  y 

should  then  increase  of  decrease  very  rapidly  as  ~|~p  increases  or 

decreases,  and  should  be  below  T*  when  is  low  and  above  T« 

when  is  high. 

The  suddeness  of  the  increase  in  Tg  as  reaches  273  K  was 

not  anticipated,  and  it  suggests  that  there  may  actually  be  two  stable 
equilibria,  with  an  unstable  euilibrium  in  between.  This  idea  is 
supported  by  the  observation  that  for  "T"<^  —  273  K  ,  the  rate  at  which 
Tif  departs  from  its  initial  value  continually  increases  for  the  first 
three  years.  To  test  our  hypothesis  we  have  performed  some  additional 
2-year  runs,  increasing  in  steps  from  258  K  to  288  K,  and  then 

decreasing  in  steps  back  to  258  K.  In  each  run  after  the  first  one 

the  initial  conditions  are  the  final  conditions  of  the  previous  run. 

We  find  indeed  that  between  -  268  K  and  "T^  “  278  K  there  are 
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multiple  equilibria,  comprising  a  cold  regime  which  is  approached  when 
the  initial  temperatures  are  low  and  a  warm  one  which  is  approached  when 
they  are  high.  Fig.  1  shows  the  equilibrium  values  of  T<*  plotted 
against  Te?  .  The  large  dots  indicate  the  results  of  particular 
runs.  The  dots  fit  two  smooth  curves. 

Presumably  these  curves  are  portions  of  a  single  smooth  S-shaped 
curve;  the  intermediate  portion  passes  through  the  unstable  equilibria, 
and  in  Fig.  1  it  has  been  sketched  by  eye.  It  could  be  located  by  a 
successive-approximation  procedure  beginning  with  a  good  initial  guess, 
such  as  the  guess  shown  in  Fig.  1. 

We  cannot  say  at  this  point  whether  the  model  has  passed  or  flunked 
its  first  test.  The  cold  and  warm  regimes  seem  unreasonably  cold  and  warm 
by  comparison  with  the  real  atmosphere,  and  indeed  the  oceans  would  be 
thoroughly  frozen  in  the  oold  regime  while  the  implicit  assumption 
that  the  mixing  ratios  are  fairly  small  would  not  hold  in  the  warm  regime. 
However,  in  reality  no  column  of  air  with  its  underlying  ocean  is  iso¬ 
lated,  since  there  is  always  a  circulation  in  the  atmosphere  and  in  the 
ocean  between  the  more  strongly  and  less  strongly  heated  regions.  Also, 
we  do  not  know  how  the  real  atmosphere  would  behave  if  the  earth's 
surface  were  all  ocean. 

The  absence  of  a  stable  equilibrium  value  of  between  260  K 

and  290  K  is  not  necessarily  serious.  The  approach  to  a  cold  or  a  warm 
stable  equilibrium  is  so  slow  that  when  Tq  varies  horizontally ,  the 
resulting  circulations  will  have  ample  time  to  prevent  either  equilibrium 
from  being  approached.  We  might  also  observe  that  a  curve  of  Tg 
against  5b  ^  would  have  a  nearly  constant  slope,  so  that  if  the 
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sea-surface  temperature  had  been  prespeci f ied,  this  phenomenon  would 
not  have  occurred. 

7.  Concluding  remarks 

We  have  constructed  a  low-order  moist  general  circulation  model 
which  contains  20  prognostic  ordinary  differential  equations  if  the 
ocean  is  assumed  to  have  zero  heat  capacity,  and  27  equations  if  it  has 
finite  heat  capacity.  The  model  atmosphere  is  composed  of  dry  air, 
water  vapor,  and  liquid  water,  anti  the  thermodynamic  and  radiative  effects 
of  water  vapor  and  clouds  are  included. 

We  have  run  the  model  for  the  special  case  when  the  solar  heating  is 
horizontally  uniform,  and  have  obtained  an  unexpected  result,  namely, 
that  for  certain  intensities  of  solar  heating  there  are  two  decidedly 
different  equilibrium  solutions.  We  have  not  yet  performed  any  runs 
with  horizontally  variable  solar  heating,  which  should  prove  to  be  of 
greater  interest. 

We  have  constructed  the  model  in  such  a  way  that  modifications  may 
easily  be  made,  and  we  anticipate  that, in  subsequent  development  of  the 
model,  modifications  will  be  made.  These  might  include  changes  in  the 
numerical  values  of  constants,  alterations  <t  specific  formulas  such  as 
the  one  relating  cloud  amount  to  relative  humidity,  or  explicit  recog¬ 
nition  of  the  partially  absorbing  wuter-vapoi  Lunds.  It  is  easy  to  add 
to  the  model,  say  by  giving  the  ocean  sut face  a  prespo<  i f led  albedo, 
including  the  absorption  and  emission  of  long-wave  radiation  by  a 
prespecified  concentration  of  carbon  dioxide,  or  including  some'  absorption 
of  short-wave  radiation  by  the  atmosphere. 
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We  must  recognize,-  however  that  the  low  horizontal  and  vertical 
resolution  preclude  from  the  start  the  possibility  of  closely  reproducing 
the  atmospheric  circulation.  The  purpose  of  the  model  is  to  test  quali¬ 
tatively  the  influence  of  various  moist  processes  over  a  far  wider 
range  of  conditions  than  is  economically  feasible  with  a  large  model. 
Inclusion  of  such  processes  as  reflection  by  the  ocean  is  desirable 
only  if  it  contributes  to  the  main  purpose. 

We  believe  that  the  weakest  features  of  the  model  are  the  assumptions 
of  a  uniform  lapse  rate  of  temperature  and  a  uniform  tropopause  pressure. 
Together  these  imply  a  perfect  relation  between  surface  temperature 
and  tropopause  temperature,  and  hence,  except  when  very  little  moisture 
and  cloudiness  is  present,  a  close  relation  between  surface  temperature 
and  the  temperature  at  which  the  bulk  of  the  outgoing  radiation  takes 
place.  There  should  therefore  be  large  variations  of  outgoing  radiation 
with  latitude.  In  t hr-  real  atmosphere  these  variations  are  not  so  great, 
because  in  the  tropics  the  uppermost  layers  of  water  vapor  are  very  high, 
and  arc  as  cold  as  those  in  the  polar  regions,  which  are  not  so  high. 

Hence  the  real  atmosphere  should  require  a  greater  cross-latitude  heat 
transport  than  the  model.  Also,  the  stabilization  of  the  lapse  rate 
when  kinetic  encrqy  is  released  would  be  a  desirable  feature. 

It  is  our  ultimate  intention  to  attempt  to  eliminate  these  weaknesses 
We  should  be  able  to  allow  the  vertically  uniform  lapse  rate  to  vary 
horizontally,  thereby  adding  seven  more;  prognostic  equations  to  the  model. 
We  prefer  to  find  some  means  for  representing  the  tropopause  height  in 
term;;  of  the  surface  temperature  and  the  lapse  rate,  and  perhaps  some 
other  quantities,  instead  of  letting  the  tropopause  height  become  still 
another  prognostic  variable. 
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